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The plasma wake field (PWF) excitation^ is one of the efficient mechanisms to generate, in a plasma, ultra-intense 
electric and magnetic fields. To produce these fields, a relativistic electron/positron beam is launched into the 
plasma. Due to the very strong nonlinear and collective beam-plasma interaction, the beam becomes the driver of a 
large amplitude plasma wave that follows the beam with almost its own speed (plasma wake) and carries out both 
longitudinal and transverse electric fields (plasma wake fields). If the beam length is much greater than the plasma 
wavelength (long beam limit), the entire beam experiences the effects of the wake fields that itself has produced (self 
interaction) . 

We study the self interaction of an electron/positron beam travelling in a magnetoactive, collisionless, cold plasma in 
the overdense regime, i.e., no 3> (no and nt being the unperturbed plasma and beam densities, respectively). We assume 
that: a strong constant and uniform external magnetic field acts along the z-axis, Bo = Boe z ; the electron/positron beam is 
travelling along the z-axis, i.e. vj, = /3ce z (j3 ~ 1); the ions are at rest to form a background of positive charge. Moreover, 
we consider a fluid model of the beam-plasma system, characterized by the electron fluid velocity, u(r, t), the electron plasma 
and beam number densities, n(r,t) and pb(r,t), respectively, the electron plasma and the beam current densities, — enu and 
f3qpt,ce z , respectively, where q — — e and q = e for electrons and positrons, respectively. We ignore the longitudinal beam 
dynamics and concentrate our investigation on the transverse effects only. We express the total electric and magnetic fields of 
the beam-plasma system in terms of the four-potential (A, cj>) and introduce small deviations of all the quantities with respect 
to the initial state n = no, u = 0, A = Ao(r) = (Bo x r) /2, = 0; then we write u = e z u\ z + uix, A = Ao + e z A± z + Aij_, 
4> — cj>i. Furthermore, we assume the long beam limit, i.e. d 2 /d£ 2 <C uj 2 /c 2 = k 2 (oj p being the electron plasma frequency), 
which under suitable boundary conditions implies that |Aix| -C A\ z and |u.ix| <C u\ z . Consequently, after linearizing the 
Lorentz-Maxwell system of equations, we obtain the equation for the dimensionless wake potential U w , driven by the beam 
density, i.e., (Vj_ — ft 2 ) U w = K 2 /9;,/no7o, where ujuh is the upper hybrid frequency, re = ui 2 /ujuhc 2 , Vi is the transverse part 
of the gradient operator, and U w = (A\ z — 4>i) /niojoc 2 with A\ z — Ai z (r±,£) and 4>i = 4>i( r ±i£.) the perturbations of both 
the longitudinal component of the vector potential and the electric potential, respectively. Here, ri is the transverse position 
vector, £ = z — pet ~ z — ct plays the role of time-like variable, mo and 70 are the rest mass and the unperturbed relativistic 
gamma factor of the single particle in the electron/positron beam, respectively. To provide a self consistent description of the 
transverse electron/positron beam dynamics, we have to consider the total force acting on the single particle of the beam, that 
accounts for the transverse gradient of U w as well as the effects produced by the exte rnal magnetic field Bo. To this end, we 
describe the transverse beam dynamics by means of the Thermal Wave Model (TWM)®. 

TWM provides an effective description of the charged-particle beam dynamics in terms of a complex function, say ty, called 
beam wave function (BWF). The transverse spatio-temporal evolution of the BWF is given by ied^ /d£ = H (rx, — ieVx, £) $ ', 
where, for normalized p&(rx,{) = (A r /a 2 )|^(rx, £)| 2 (N and a z being the total number of particles and the beam length, 
respectively)^, H(r±, p ± , £) is the effective Hamiltonian describing the perturbed transverse motion of a single particle of the 
beam, and e the transverse beam emittance. Provided that the longitudinal dynamics is ignored, i.e. p = e z mo7oc + p 1± , it is 
easily seen that, "H = AH/Ho = (H—Ho)/Ho is the relative first-order perturbation of the single electron/positron Hamiltonian 
H(r,p,0 = cy/(p- g -A) 2 + m 2 c 2 + q( />, where H = mo7oc 2 is the initial unperturbed total energy of the single particle of the 
beam. Consequently, the above equations for BWF and U w can be cast as the following Zakharov system of equations, viz., 
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where we have put $(rx, f, £) ~ exp [im (<p — fc c £/2)] ip m (r±, f) with m integer and K — (uj c /2 r yoc) 2 = (qBo/2mo r yoc 2 ) 2 = 
(k c /2) 2 . This system governs the self consistent spatio temporal evolution of the PWF self-interaction of the electron/positron 
beam. In principle, once eq.(p| is solved for U w , we get the functional U = U [iV'ml 2 ] that makes eq. Jill the generalized 
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nonlinear Schrodinger equation (NLSE). Note that, due to the helicity phase factor exp [im (ip — fc c £/2)] for m 7^ the BWF 
describes vortices states associated with the orbital angular momentum of the beam particles in the external magnetic field 
(m plays the role of vortex charge). Based on the pair of eqs. |T| and an investigation, both analytical and numerical, 
has been carried out, taking into account the diverse limiting cases. Hereafter, and in RefN, we give a summary of the maim 
results. 

First of all, in the linear limit, i.e. U = U \ \^>m\ 2 — (the beam is travelling along Bo and the self interaction is assumed 
negligible), for arbitrary integer m, the pair of eqs. (I]) and |2| reduces to a linear Schrodinger equation that admits solutions 
in the form of Laguerre-Gauss modes, say tp p ,m(f±, £ji (p is an arbitrary integer), whose transverse size, cr p ,m,(£), * s proportional 
to the one associated with the purely Gaussian mode (m = p = 0), say cr(£). The latter satisfies the following envelope 
equation: d 2 o/d£ 2 + Ka — e 2 /er 3 = 0. For a long enough region along z, it describes the typical sausage-like transverse beam 
size modulations (betatron oscillations), preserving the collapse (focusing to a single point). Fig. [I] shows the density plot for 
|^ P ,m| 2 in the normalized plane (x/cr, y/o-) for some combinations of p and m. For a given m, the number of rings increases 
as p increases. When m = 0, the structure of the density plots is always constituted by a central core plus p rings, whilst for 
ra/0, the density plots are associated with the diverse states of vortices. 
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FIG. 1: Left: Density plots for |\& m ,p\ in the normalized plane (x /a , y /a) for some combinations of p and m. Right: qualitative 
representation of the BWF phase for different values of the vortex charge m when p — 0. 



When the self interaction is not negligible, one can suitably investigate two limiting cases, 
(i). If the transverse beam size is much greater than the plasma wavelength, i.e. k v a 3> 1, the system ([T]) and Q reduces to 
the 2D Gross-Pitaevskii equation, i.e., 
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In RefP we carry out a numerical analysis of this sort of NLSE, showing the existence of vortices, nonlinear coherent states 
(2D solitons) and beam halos. here, from the virial equation associated to this NLSE, we find the following envelope equation: 

d 2 a 2 m /d^ 2 + 4Ka 2 m = 4A m , 

where A m = 7re 2 [\d^ m /dr ± \ 2 + (m 2 /r\) \i> m \' 2 ] r ± dr ± -(nN/a z n 'y ) J °° | ip m | 4 r ± dr ± + (K/2)a 2 n is constant of motion, 
i.e. dAm/d^ = 0, and a m = [27r r\ \ipm{r±, £,)\ 2 r± dr±j 1 2 is the transverse beam size. Note that the interplay between 
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the positive and the negative terms in the definition of Am can make this quantity positive, negative or zero. If Am > 0, in 
correspondence of the initial conditions ao = cr m (0), and a'o = (dam/d^)^^ , the envelope equation describes stable oscillations 
in a m with frequency 2vA, in the range < |oq 2 + \Ko~ 2 - ) < Am, whilst a m would reach zero in a finite time in the range 
< Am < |o"o 2 + \Koq. If Am < 0, always a m would reach zero in a finite time. The latter cases would correspond to a 
collapse instability. However, as a m reaches very small values, the condition k p a m 3> 1 is no longer satisfied and the collapse 
does not take place. It is interesting to describe the self-interaction of an electron/positron beam that enters a thin plasma 
slab [plasma lens) of length Z at £ = where, as initial conditions, is assumed that a' = and the BWF is purely Gaussian, 
ip(rx,0) = exp[— r 2 L /2<T 2 ]/y / 7Fao. Provided thatx/AZ <C 1, the self-interaction is very short, namely the the BWF remains 
almost unchanged except for the appearance of a chirping phase factor, that accounts for a strong change in the particle 
momentum distributions in the transverse plane (kick approximation), viz. ,ip(r±, I) ~ ip(r±, 0) exp[ir 2 _/2ep(l) + icj>(l)], where 
p(£) and <j)(0 are the curvature radius of the wavefront and a homogeneous phase at location £, respectively, such that, as 
£ — > 0, p — > oo and <f> — > 0. At the lens exit (£ = I), the transverse beam size is a (I) ~ \J a \ + 2 (Ao — Kgq) Z 2 (Aq being 
the constant of motion corresponding to the Gaussian beam). Consequently, under the condition Ao — Kcrfi < 0, the beam is 
focussed and out of the lens it reaches a minimum spot size that is greater than zero, according to the envelope equation in the 
vacuum (i.e., no plasma and K — 0), viz., d 2 a/d^ 2 — e 2 /<r 3 = 0. We have weak focusing (strong focusing) if the above condition 
is satisfied for positive (negative) values of Ao. 

(ii). If Zcp<7 <g 1, eqs. |TJ and (J2J reduces to a nonlocal NLS equation whose aberrationless (i.e. Gaussian) approximate 
solution leads to the envelope equation d 2 a/d£ 2 + Ka + n/a — e 2 /a 3 — 0, where 

N/mo'yo^uH '^)- It admits the 

self- equilibrium solution a' eq = \J~n + y/l + (J^ o~ eq , where a eq = e 1//2 /A 1//4 is the corresponding self-equilibrium solution for 

the linear case n — and fj, = rj/2e\ /r K. It is easily seen that a' eq < a eq in all the range < pi < oo. This implies that the term 
proportional to 1/a accounts for the squeezing of the beam. 
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